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Abstract 


The corrosion inhibition performance of Hexa (3-methoxy propan-1,2 diol) 
cyclotriphosphazene (HMC) on carbon steel in 3% NaCl solution was 
investigated by weight loss (WL), potentiodynamic polarization (PDP), 
electrochemical impedance spectroscopy (EIS) measurements, Density functional 
theory (DFT) and Monte Carlo (MC) simulation. The corrosion inhibition 
efficiency at optimum concentration (107°M) is 99% of HMC at 298 K. The 
corrosion inhibition efficiency at 107° M decreases with increase in temperature. 
The adsorption of HMC on the surface of carbon steel obeyed Langmuir 
isotherm. Potentiodynamic polarization study confirmed that inhibitor anodic- 
type. DFT and Monte Carlo (MC) simulations based computational approaches 
were under taken to support the experimental findings. DFT studies revealed that 
HMC interact with metallic surface through donor-acceptor interactions in which 
the anionic parts act as electron donor (HOMO) and cationic parts behaved as 
electron acceptor (LUMO). The MC simulations study showed that studied HMC 


adsorb spontaneously on Fe (110) surface. 


Keywords: Materials chemistry, Theoretical chemistry, Inorganic chemistry 


1. Introduction 


Steel alloys are widely utilized in many applications including petroleum industries, 
especially carbon steel are prone to corrosion when they come into contact with cor- 
rosive materials, especially in an environment containing Cl ions. Presence of CI- 
ion on the metallic surface initiates the formation of metal oxides; thereby adversely 
affect the physiochemical properties of the metals. The use of corrosion inhibitors, 
mainly heterocyclic compounds have been reported to be one of the most effective 
methods for against corrosion because they can easily adsorb on metallic surface via 
their molecular structures P, O, S, N, 7t-electrons and aromatic rings [1, 2, 3, 4, 5]. 
Many of these heterocyclic compounds are obtained efficiently using easy-access 
synthesis pathways and are usually used at lower concentrations during the corrosion 
inhibition process [6, 7, 8]. Among this group of heterocyclic compounds, we find 
Hexachlorocyclotriphosphazene. Hexachlorocyclotriphosphazene are an important 
group containing the formula (NPC1,)3. The molecule has a cyclic backbone consist- 
ing of alternating P and N atoms. They have active phosphorus-halogen bonds that 
can be replaced with nucleophiles, leading to the formation of different types of 
product [9, 10, 11]. Cyclotriphosphazene is well established nitrogen rich ring con- 
taining molecule with three additional nitrogen atoms those can easily protonate and 
enhance the solubility of polar solvents. Recently, derivatization of cyclotriphospha- 
zene gained the significant advanced for variety of purposes including corrosion in- 


hibition. Cyclotriphosphazene have been investigated as effective corrosion 
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inhibitors for metals and alloys in aggressive solutions owing to their high protection 
ability which is in turn attributed to the adsorption of these compounds by their non- 
bonding electrons of nitrogen atoms and 7t-electrons of three double (-P=N-) bonds. 
In this study demonstrated the effect of cyclotriphosphazene derivative designated as 
Hexa (3-methoxy propan-1,2-diol) cyclotriphosphazene (HMC) having either six 
identical oxygen containing 3-methoxy propan-1,2 diols. Last decade has witnessed 
a dramatic increase of using Density Functional Theory methods (DFT) to magnif- 
icently define the chemical reactivity of inhibitors and their adsorption competence 
on metal surface [12, 13, 14, 15, 16, 17, 18]. Such computations have been widely 
used to analyze the molecular electronic structures of adsorption-type inhibitors us- 
ing a number of quantum chemical descriptors which gives important insights on 
corrosion inhibition mechanisms. Coupled with the fact that metal inhibitor interac- 
tions can be best describe with the aid of the Hard and Soft Acid and Bases (HSAB) 
principle [19], where organic inhibitors can be regarded as Lewis bases that offer un- 
shared non-bonding and 7t-electrons to the empty d-orbitals of the iron atoms that 
acted as electron acceptor (Lewis acid). Therefore, inhibition behavior of the organic 
inhibitors including HMC is largely influenced by the electronic properties of the in- 
hibitor molecules. Monte Carlo (MC) simulation has been another effective tool to 


study the interaction of inhibitors with the metal surface or similar problems. 


In view of this, present investigation deals with the implementation of one cyclotri- 
phosphazene derivatives synthesized by condensation of Hexachlorocyclotriphos- 
phazene and 2,3-epoxy-1-propanol and demonstrated as corrosion inhibitor for 
marine-environment (3% NaCl) corrosion of carbon steel. Several commonly em- 
ployed methods such as chemical (weight loss), electrochemical impedance spec- 
troscopy (EIS), and potentiodynamic polarization (PDP) have been undertaken to 
study the inhibition ability of the Hexa (3-methoxy propan-1,2-diol) cyclotriphos- 
phazene (HMC). The adsorption and inhibition behavior of HMC were supported 
by density functional theory (DFT) based quantum chemical calculations method. 
A significant correlation among the parameters of experimental and theoretical 
methods has been achieved in present study. The most stable configuration and 
adsorption energies for the interaction of the inhibitor with Fe (110) surface were 


studied by Monte Carlo simulation. 


2. Experimental 
2.1. Materials 


Carbon steel specimens were used as working (test) material. The corrosive medium 
of 3% NaCl concentration was manufactured by dilution with distilled water. The 
inhibition effect of HMC was tested at four different concentrations denoted by, 
107°, 107, 1075 and 1076 M. The inhibitor, HMC was synthesized according a re- 


ported procedure [6]. In brief, a mixture of an (2,3-epoxy-1-propanol) (2) compound 
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(2g) and toluene (14.25mL) containing 2.22 g of triéthyl amine at room temperature. 
The mixture was cooled (ice-water). After was added 1.5 g of hexachloro cyclotri- 
phosphazene (1) and toluene (7.5 mL) in two-necked round bottom flask, fitted with 
a condenser and a magnet stir bar was added drop wise for 3 h. After the addition was 
completed, the stirring was continued again for another 45 h at room temperature. 
The residual solvent was removed using a rotary evaporator in order to obtain the 
desired brown viscous solution of resin (3). After we take 5 g of product obtained 
epoxy resin (ER) (3) and added 10 mL of basic solution containing 2.1 g of 
NaOH. The mixture was neutralized, filtered and dried by NazSO, and to remove 


water. 
Schematic diagram for the preparation of the HMC is shown in Fig. 1. 


The structure of the produced was verified by 'H NMR and FT-IR spectroscopic 
characterization techniques. HMC: brown solid product, very soluble in water, yield 
95%; 'H-NMR (CDCI, 300 MHz): dppm = 2.7 (dd, 2H, -O-CH2), 3.38 (m, 1H, 
CH), 4 (dd, 2H, CH).7'P-NMR: ôppm = 9.32 (s). FTIR (cm~'):1266 and 1200 
(P=N stretching), 1013 (P-O-C), 2927, 2965 (asymmetric and symmetric y of 
C—H units), 2800, 3600 (—OH). 


2.2. Methods 
2.2.1. Weight loss measurements 


The weight loss analysis was carried out using the carbon steel of previously re- 
ported composition and procedure as described in our earlier studies [20, 21]. The 


gravimetric parameters were derived using following equations [22]: 
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Fig. 1. A representative scheme showing the total synthesis of Hexa (3-methoxy propan-1,2-diol) cyclo- 
triphosphazene (HMC). 
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corr = — 1 
W. F (1) 
-mw 5 
Wo 
-w, 
nwi, = ——— x 100 (3) 


0 


where Wo and W; represent the weight loss values without and with HMC, respec- 


tively, t(h) is the immersion time and S is the total surface area of the metal sample. 


2.2.2. Electrochemical measurements 


Instruments and experimental procedures were similar as described in our previous 
reported literatures [23, 24, 25]. 


From EIS and polarization techniques, inhibition effect of HMC was evaluated using 
Eqs. (4) and (5). 


The inhibition efficiency ngrs% is calculated from resistance polarization R, of HMC 


given by Eq. (5): 


Rp — R? 
Mers (%) =—=— x 100 (4) 
Rp 
leoi 
Nppp % = (1 -y ) x 100 (5) 


where, R? and R, represent the resistance polarization values without and with 
HMC, respectively. Whereas, SA and isor represent the corrosion current density 


values without and with the inhibitor respectively. 


2.2.3. Computational details 


DFT on the HMC in its neutral as well as solvated phase was carried out using 
Gaussian 09 similar to our previously reported procedures [26, 27, 28, 29, 30, 31, 
32,33]. 


Several DFT parameters were computed using following relationships [34, 35, 36, 
37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48]: As per the theorem of Koopman, 
the values of frontier molecular orbitals can be correlated with ionization potential 
(I) as well as electron affinity (A) as per the Eqs. (6) and (7). The energies the en- 
ergies are highest occupied molecular orbital (Enomo) and lowest unoccupied mo- 
lecular orbital (Eumo) and energy gap (AE = FE, umo-Exomo) have been calculated. 
Employing the magnitude of A and I, several other DFT indices like electronega- 


tivity (x), global harness (ņ) and softness (7) were computed for the forms of 
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HMC as per the Eqs. (8)— (10). It is well documented that metal inhibitor bindings 
involve donor-acceptor phenomenon. The fraction of electron transfer (AN) by 
HMC molecules was computed as per the Eq. (11). Which is the initial molecule- 
metal interaction energy (Ay) and the inhibitor molecule and the metal surface as 
per the Eqs. (12) and (13). 


I = —Exomo; (6) 
A= —E\vmo; (7) 
I+A 
x= (8) 
I-A 
1 
o=- 10 
. (10) 
Pre — Xinh 
AN = OO 11 
2(Nre + Ninn) ( ) 
(Xre ~~ Xin) 
AY = ay 12 
2(Nre + Ninh) ( ) 
AE pack-donation = -1 (13) 


Where, where ¢ corresponds to the work function, which is used as the appropriate 
measure of electronegativity of iron, and np. = 0. The theoretical value of > is 4.82 
eV for Fe (110) plan which is reported to have the higher stabilization energy. 
Eq. 13implies that when n > 0 or AE,.a< 0, back-donation from the molecule to 


metal is energetically favored. 


It is very necessary to identify the reactive centres present within the inhibitor mol- 
ecules. The Fukui Indices (FIs) analyses help to determine and recognize the sites of 
the molecule which are responsible for nucleophilic or electrophilic reactions. The 
presence of the nucleophilic or electrophilic center facilitates the interaction of the 


inhibitor molecules with the surface atoms of the metal. 


Fukui functions for nucleophilic and electrophilic attacks are obtained from finite 


difference approximations as the follows: 


JEC) =p.(N+1)—p,(N) For nucleophilic attack (14) 
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f(r) =p,(N) —p,(N—-1) For electrophilic attack (15) 


px(N),px(N-1) and p(N+1) are the gross electronic populations of the site k in 


neutral, cation and anion, respectively. 


Monte Carlo simulations help us to find the most stable adsorption sites on metal and 
metal oxide surfaces through finding the low energy adsorption sites on both peri- 
odic and non-periodic substrates or to investigate the preferential adsorption of mix- 
tures of adsorbate components. In order to investigate the most suitable adsorption 
modes for the neutral inhibitor, the HMC molecule and Fe20; surface were prepared. 
In this study, the HMC molecule was drawn and the hydrogens were adjusted, then 
the molecules were cleaned using sketch tools available in the Materials Visualizer. 
The sketched and fully optimized geometry structure of HMC is shown in Fig. 8. 
The molecule that are shown in Fig. 8 were optimized (i.e., energy minimizing) in 
the Monte Carlo simulations using the forcite module and the Condensed- Phase 
Optimized Molecular Potentials for Atomistic Simulation Studies (COMPASS) 
force field. For cite is an advanced classical mechanics tool that allows energy cal- 
culations and geometry optimizations. COMPASS force field which used to opti- 
mize the structures the HMC of all components of the system of interest (Fe* 
inhibitor) and represents a technology break-through in force field method. COM- 
PASS is the first ab initio force field that enables accurate and simultaneous predic- 
tion of chemical properties (structural, conformational, vibrational, etc.) and 
condensed-phase properties (equation of state, cohesive energies, etc.) for a broad 
range of chemical systems. It is also the first high quality force field to consolidate 
parameters of organic and inorganic materials. Materials Studio provides the unit 
cell structures of metal oxides with the associated experimental lattice parameters. 
Metal oxide surfaces were prepared by employing the desired cleavage planes hkl 
(110), Using the surface builder module of Materials Studio. The Metal oxide sur- 
face must be large enough to accommodate the inhibiter molecules; therefore, the 
constructed surface was enlarged by periodic replication to triple U xV (3 x 3) 
in order to expose a more realistic surface area for docking the inhibitor molecules. 
The surface with parameters of 45.7 and 27.9 A was produced. It is important that the 
size of the vacuum is great enough such that the non-bonded calculation for the 
adsorbate does not interact with the periodic image of the bottom layer of atoms 
in the surface; thus, a vacuum slab with thickness 31.7 A was built using Crystal 
Builder. A low-energy adsorption site is identified by carrying out a Monte Carlo 
search. This process is repeated to identify further local energy minima. During 
the course of the simulation, adsorbate molecule are randomly rotated and translated 
around the substrate. The configuration that results from one of these steps is 
accepted or rejected according to the selection rules of the Metropolis Monte Carlo 
method. The Monte Carlo (MC) simulation was performed using Materials Studio 
5.5 program [49, 50, 51, 52, 53]. 
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3. Results and discussion 
3.1. Weight loss measurements (WL) 


The WL parameters derived for steel corrosion after 24 h immersion time in 3 % 
NaCl solution is listed in Table 1. Visualization of the results showed that inhibition 
performance of HMC is concentration dependent it increases on increasing the con- 
centration of HMC and maximum efficiency was observed at 107°? M. The increase 
in the inhibition effectiveness of HMC with their concentration is might be attributed 


to the increase in surface coverage values at higher concentrations. 


The process of corrosion inhibition was highly influenced by the temperature. For 
the purpose of further understanding the effect of temperatures, the corrosion inhibi- 
tion ability of HMC in 3% NaCl solution was studied at various temperatures in the 
rang 298—328 K for 24 h of immersion. The evolution of W.o, and nw, with 
various temperature for carbon steel in 3% NaCl solution in the without and with 
addition 107° M of HMC is given in Table 2. 


The results indicate that the values of nw, obtained from the weight loss for carbon 
steel corrosion in 3% NaCl solution in the without and with addition 107° M of 
HMC at different temperatures decreased with increasing temperature. The 
maximum value of nw ,% obtained for 107° M of HMC is 96 % at 298 K. This 


behavior indicates desorption of inhibitor molecule. 


3.2. Potentiodynamic polarization measurements (PDP) 


The PDP curves for carbon steel in 3 % NaCl are shown in Fig. 2. Parameters are 
presented in Table 3. Carbon steel readily undergoes anodic metallic dissolution 


and cathodic hydrogen evolution reactions using following schemes. 


Anodic dissolution reactions [54]: 


Fe Fe + 2e- (16) 


Table 1. W..,,, and nw obtained from weight loss measurements for carbon steel 
in 3% NaCl with and without various concentration of HMC at 298 K for 24 h. 


Inh cC M) Weorr (mg/em*h’) nwi(%) 
Blank = 3 = 
10°° 0.23 92 
HMC 10 0.22 93 
1074 0.17 94 
107° 0.12 96 
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Table 2. Effect of temperature on the corrosion rate (W,,,,) and corrosion current 
densities (igo) for carbon steel in 3% NaCl solution in the absence and 
additionl10~* M of HMC. 


T (K) Weorr (mg/cm?h') İcorr nwi (%) nppr (%) 
uA/Cm? 


Blank 107° M of HMC Blank 107° M of HMC — — 


298 3 0.12 280 3 96 99 
308 5.51 0.41 345 23.4 92 93 
318 8.89 1.31 556 85.8 85 85 
328 11.45 2.05 706 137.5 82 80 


Fe +CT © (Fecr) 


ads 


(FeCl)... (FeCl) as + le (18) 
(FeCl) a FeCl’ + 1e (19) 
FeCl @ Fe** + Cr (20) 


Cathodic hydrogen evolution reactions: 
O2 +2H20 + 4e > 40H- (21) 
2H,O + 2e —H,+20H- (22) 


However, in its presence HMC forms protective inhibitive barrier which separates 


metal from the aggressive 3 % NaCl solution and protect from corrosion. 


i (A/cm?) 


=== 10* M of HMC 
—=< 10° M of HMC 
42 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 02 
E (W/Ag/AgCl) 


Fig. 2. PDP curves for carbon steel corrosion in 3% NaCl solution containing different concentrations of 
HMC. 
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Table 3. PDP parameters for carbon steel corrosion in 3% NaCl solution with and 


without addition of different concentrations of HMC. 


Inh CM) Eor (MV/Ag/AQCl)  icorr (HA/em”) Ba (VTD -6e (V~") neor% 


Blank — —559 280 0.12 0.12 — 
1076 —552 40 0.26 0.19 85 
HMC 10 —498 27 0.22 0.18 90 
1074 —463 14 0.29 0.16 94 
107° —320 3 0.26 0.03 99 


From the results obtained, the consequences is was perceived that icorr on the two 
cathodic and anodic polarization branches decrease in the presence of the HMC. 
Correspondingly, the inhibition efficacy increases by increasing the concentration 
of the HMC. The maximum efficiency of 99 % was obtained at 107°? M. Further- 
more, as shown in Fig. 2 and Table 3, with the addition of HMC molecules shifts 
corrosion potential is shifted to more anodic potential and also shows a significant 
change the anodic polarization branches. Anodic dominance of the HMC can be 
easily observed from the results presented in the table [55]. These results suggest 
that HMC molecules reduces the anodic metallic dissolution and retards the reduc- 


tion of oxygen. 


The influence of temperature on the various corrosion parameters icor and Nppp% 
obtained from the PDP of carbon steel was studied in 3% NaCl solution in the tem- 
perature rang (298—328 K)without and with addition 107° M of HMC. Results are 
listed in Table 2. 


From the results of Table 2 shown that, as the temperature increased in in 3% NaCl 
solution in the without and with addition 107° M of HMC, the values of icorr increase 
and Nppp% decrease. This result can be explaining by desorption of HMC molecules 


from carbon steel surface. 


3.3. Electrochemical impedance spectroscopy measurements 


Nyquist plots for the carbon steel electrode in 3% NaCl solution and containing 
various concentrations of HMC after 30 min of immersion time are presented in 
Fig. 3. 


It can be seen that in the absence of inhibitor, Nyquist plot contain only one semi- 
circle. This indicates that corrosion on the carbon steel surface is mainly controlled 
by the charge transfer process [56, 57]. In Nyquist plot, difference in the lower fre- 
quency and higher frequency in the real axis is generally considered as charge trans- 
fer resistance (Re). The charge transfer resistance must be consistent with the 


resistance between metal and Outer Helmholtz Plane (OHP). Thus, the resistance 
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= Blank 
> 10° MofHMC 


-Z" (Q.cm’) 


Z (Q.cm’) 


Fig. 3. Nyquist diagrams for carbon steel in 3% NaCl solution containing different concentrations of 
HMC. 


coming from the metal/solution interface in absence of the inhibitor is mainly due to 
the charge transfer resistance (R.), diffuse layer resistance (Ra) and accumulation 
resistance (R,). Therefore, along with Res Ra and R, have to be taken in account. 
Thus, difference in the lower frequency and higher frequency in the real axis is 
considered as polarization resistance (Rp). It is also seen from the Nyquist plot 
that the diameter of semi-circle is concomitantly increased in the presence of the in- 
hibitor. This phenomenon suggests the formation of protective layer on the metal 
surface. Thus, in presence of inhibitor film resistance (Rp) have to be considered 
in the polarization resistance (R,). So, in presence of inhibitor, R, is designated as 
= (Ra + Ra + Ra + Rp). The addition of HMC to the aggressive solution leads to 
a change of the impedance diagrams in both shape and size, which a depressed semi- 
circle at the high frequency part of the spectrum was observed and a second time 
constant appeared at low frequency region. The first capacitive loop appearing at 
high frequency region was attributed to the charge transfer resistance and the diffuse 
layer resistance, the second one at low frequencies related to the adsorption of inhib- 
itor molecules on the metal surface and all other accumulated kinds at the metal/so- 
lution interface (inhibitor molecules, corrosion products, etc.). As seen from Fig. 3, 
the R, values increased with the HMC concentration, which can be attributed to the 
formation of a protective over layer at the metal surface, which becomes a barrier for 
the mass and the charge transfers. The double layer can be represented by the elec- 
trical equivalent circuit diagrams to the model metal/solution interface. The corre- 
sponding electrical equivalent circuit model for the carbon steel in uninhibited 
solution is given in Fig. 4c with a schematic representation of the potential distribu- 
tions on the metal/solution interface (Fig. 4a) and the resistances related to the dou- 
ble layer (Fig. 4b). According to Fig. 4, the R, includes the Res, the Ra and all the 


other accumulated kinds of ions, corrosion products, etc. at metal/solution interface. 
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Metal 
Solution 


Solution 
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(b) 


1 i | 
| i 
| I 
l | 
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metal/soltion interface 


1 Rp =RettRatRa 


Fig. 4. The potential distributions on the metal/solution interface (a), the resistances related double layer 


(b), and proposed electrical equivalent circuit diagram for blank solution (c). 


The Ra must be corresponded to the resistance between the metal and OHP (outer 
Helmholtz plane). All the other resistances must be considered as the simple resis- 
tance on the current crossing away. In the case of the inhibited solutions, related elec- 
trical equivalent circuit model is given in Fig. 5b with a schematic representation of 
metal/solution interface (Fig. 5a) [58, 59, 60]. 


The obtained Nyquist plots have been fitted using these equivalents circuits where 
the constant phase element and the polarization resistance are in parallel combination 
and these two are connected in series with the electrolyte resistance (R,). The con- 


stant phase element is defined by the following Eq. (23) [59]: 


1 


Olay (23) 


Zcre = 


Where, Q is the proportionality coefficient, w represents angular frequency, i is the 
imaginary unit and n is a measure of surface irregularity. For whole number of n = 1, 


0, -1, CPE represent a capacitance (C), resistance (R) and inductance, respectively. 


(a) 


Pore or free surface 


Metal 


Rp =Rf+Rpor RporRet+Ra+Ra 


Ct : film capacitance Rr: film resistance 


Fig. 5. Schematic representation of metal/inhibitor/solution interface (a) and proposed equivalent elec- 


trical circuit diagram for EIS results of the carbon steel in inhibited solutions (b). 
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A direct correlation between the polarization resistance and double layer capacitance 


(Ca at the metal/solution interface is obtained by the following Eq. (24): 


1/n 
Q 


Ca — RT (24) 


The electrochemical parameters values are listed in Table 4. Results depicted in 
Table 4 shows that R, value increases with increasing the concentration of HMC. 
It may suggest that a protective layer has been formed on the electrode surface. 
This layer creates a barrier for charge and mass transfer at the metal-solution inter- 
face [59]. It is also observed that with increasing inhibitor concentration the Car 
values are decreased, which can be described by an increase in the thickness of elec- 


tric double layer or by the decrease of local dielectric constant [59]. 


3.4. Adsorption studies 


The adsorption isotherm is frequently utilized to inspect the adsorption mechanism 
between the inhibitors and metal atoms at the metal surface. The most utilized 
adsorption isotherm is Langmuir [61]. The Langmuir adsorption isotherm is shown 
in Eq. (25). 


= — + Cinh (25) 


where @ is the surface coverage degree and Kaas is the adsorption process. 


Langmuir isotherm and the experimental data have a good linear relationship for the 
adsorption behavior of the investigated HMC, representing that adsorption of HMC 
on metallic surface comply with Langmuir adsorption isotherm. The association be- 
tween Cinn and Cinn/O yielded straight lines with intercept of Kaas as shown in Fig. 6. 
From the K,q, values, the AG a, are calculated using Eq. (26) [62]: 


Table 4. EIS parameters for carbon steel corrosion in 3% NaCl solution with and 


without addition of different concentrations of HMC. 


Inh C R, Rp Ca (HF/cmô) ners% 
M) (kQ.cm?) (kQ.cm?) 

Blank = 10 0.165 — — 
1076 15 1.0 165 85 

HMC 105 12 1.6 141 90 
1074 12 2.7 68 94 
107° 8 9 52 98 
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1,2x10° 


= HMC 
1,0x10° Linear Fit of HMC 


4 


8.0x10" 
6.0x10* 


4 


4.0x10" 


C „/ © (mol/L) 


4 


2.0x10" 


0.0 


-2,0x10* 


0.0 2.0x10f 4,0x10° 6,0x10* 8.0x10* 1,0x10° 


C „ mol/L) 


Fig. 6. Langmuir adsorption isotherm on weight loss of carbon steel in 3% NaCl solution containing 
different concentrations of HMC after 24 h at 298 K. 


1 —AG. 
K, <= ads 2 
i 55.3°?( RT ) (26) 


where R is the universal gas constant (8.314 J/mol. K), value 55.55 is the water con- 
centration in the solution (mol/L). In general, a high value of Kaas and low value of 
AG as indicating a strong interaction and a strong adsorption on the metal surface 
[63]. The calculated value of Kaas (1.18 Mx 10°)) and AG, (—44.57 kJ mol`’) 
indicates that HMC exhibits a strong interactions with carbon steel surface (chem- 


isorption) [64]. 


3.5. Effect of temperature 


Generally, increase in the temperature causes significant decrease in the protection 
ability of the inhibitor molecules. The influence of temperature can be best explained 
with the aid of Arrhenius Eqs. (27) and (28) [65, 66]: 


Ea 
icon = A.exXp ( — a) (27) 


E, 
Weor = Acxp( — =) (28) 


Arrhenius plots is shown in Fig. 7. Moreover, the Arrhenius equation can be con- 
verted Eqs. (29) and (30) (transition state equations) as shown [67, 68]: 


> RT (AS AH, 
icon = Fy exp( z Jeæ( RT ) (29) 
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Fig. 7. The relationship between Ln (Wor), Ln (icorr), Ln (Weorr/T) and Ln (icor/T) as a function of 1/T in 
3% NaCl solution without and with addition 107° M of HMC. 


RT AS, AH, 
Woor = IN exp( R )exr( ma ) (30) 


A is the Arrhenius exponent, T is temperature, N is Avogadro’s constant and other 
symbols have their usual meaning. Several parameters such as Ea, AH, and AS, 
calculated using Arrhenius and transition state plots (Fig. 7) are presented in 
Table 5. 


The positive value of AH, reveals the endothermic nature of metal-inhibitor interac- 
tions [69], whereas increased value of AS, suggested that replacement of the water 
molecules from ER that results into the corresponding increase in the disordering 


(randomness) on the metallic surface. 


Table 5. Values of Ea, AH, and AS, for carbon steel in 3% NaCl solution without 
and with addition 107° M of HMC. 


Method CM) E, (kJ/mol) AH, (kJ/mol) AS, (J.mol ~! K7 5 
WL Blank 35.82 39.22 —66.06 
107° M of HMC 78.95 81.52 +49.53 
PDP Blank 26.40 23.82 —118.40 
107° M of HMC 79.50 76.45 +26.25 
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3.6. Theoretical studies 
3.6.1. Quantum chemical calculation 


The anticorrosive efficiency of the Hexa (3-methoxy propan-1,2-diol) cyclotriphos- 
phazene (HMC) was investigated by DFT and MC simulations. It is worth 
mentioning that the corrosion effect of chemical inhibitors is accompanied to its 
adsorption via electronic interaction with the metal’s surface. Recently, quantum 
chemical calculation has been used to provide the scientists more information about 
structural properties of the chemical inhibitors. On the other hand the adsorption 
mechanism can be used to study the chemical reactivity of the inhibitors. The quan- 
tum chemical indicators, such as, Eyomo, ELumo, Energy gap, AE (Eyomo. 
—E,umo), electronegativity, global softness, chemical hardness and softness and 
the number of electrons transferred are very important and useful parameters to study 
the anticorrosive effect of the investigated molecule. Table 6 presents our computed 
results of the considered inhibitors in its neutral and protonated forms in both gas and 
aqueous phases. Fig. 9 shows the optimized geometrical structures and the FMO sur- 
faces (HOMO and LUMO) while the electrostatic potential maps (ESP) are pre- 
sented in Fig. 7. 


A close inspection of the optimized structure of the investigated inhibitor (Fig. 9) 
indicates that the HOMO and LUMO surfaces of the non-protonated species are 
distributed almost over the central part of the molecule except of the HOMO 


Table 6. Optimized energy (in hartree), dipole moment (11*) in Debye and 
Quantum chemical parameters in eV of ER and HMC in both gas phase and 


aqueous solution. 


Gas Solution 

ER HMC ER HMC 
Energy (a.u.) —2795.389255 —3254.345223 —2795.422351 —3254.381590 
Exomo(eV) —7.211 —6.713 —7.481 —7.370 
Eumo (eV) 0.752 0.476 0.631 0.665 
I 7.963 7.189 8.112 8.036 
A —0.752 —0.476 —0.631 —0.665 
AE 7.963 7.189 8.112 8.036 
x 3.606 3.357 3.741 3.685 
n 4.358 3.833 4.371 4.350 
o 0.229 0.261 0.229 0.230 
AN110 0.139 0.191 0.139 0.191 
Ay 0.661 0.866 0.608 0.631 
AEb-a —1.089 —0.958 —1.093 —1.088 
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structure of the non-protonated species, in which a distribution was shown over the 
one of the epoxy groups. This means that the inhibitor in its both forms has a good 
ability to donate electron to an acceptor (metal surface). The Eyjomo is a measure of a 
molecule’s tendency to receive an electron from donor species. It is well known that 
lower ELumo and/or higher Eyomo benefit(s) higher corrosion inhibition strength 
[17, 18, 70, 71]. It is clearly obvious from Table 6 that the Eyomo of the protonated 
inhibitor in aqueous medium is found to be higher than of the noon-protonated spe- 


cies. In contrast, the value of ELumo is not in agreement with our finding. 


Theoretically, another important to characterize the chemical reactivity on an inhib- 
itor is the energy gap (AE,,,). Energy gap shows the capability of inhibitor to be an 
effective corrosion inhibitor or not [72, 73]. It was found that when the AEgap de- 
creases, the adsorption efficiency between the metal surface and the inhibitors in- 
creases, because the energy to remove an electron from the HOMO must be lower 
[17, 18, 74]. 


Table 6, show that the energy break values for the non-protonated (ER) and proton- 
ated HMC species are 8.112 and 8.036 eV, respectively. The protonated form has 
inferior energy break (AE,,,) compared to the non-protonated one. This finding in- 
dicates that the adsorption performance on the steel surface of the non-protonated 


form is energetically favored. 


Another very important, both experimental and theoretical, quantum chemical 
descriptor is the chemical hardness, which describes the resistance to electron cloud 
polarization or deformation of chemical species. Additionally, chemical hardness is 
a measure of the stability of the chemical molecule, which has a high energy gap. 
On the other hand, soft molecule, which has low energy gap can be good corrosion 
inhibitor in view of the fact that soft molecule can very readily give electron of 
HOMO to metals [17, 18, 70, 71]. An inspection of Table 7 shows that the proton- 
ated species in aqueous medium is less harder and less softer than the neutral 


species. 


The value of AN reveals the capability of the examined inhibitor to transfer its elec- 
trons to metal if AN > 0 and vice versa if AN < 0. Based on this principle, it is 
apparent from the results in Table 6 that the tested inhibitor in its two forms in 
both gas and aqueous phases has higher propensity to give its electrons to a metal 


surface. 


Kovacevic and Kokalj [75] introduced another important quantum parameter, which 
is the initial molecule-metal interaction energy (Ay). In our work, Ay has been 
calculated for the investigated inhibitor in its neutral and protonated forms in both 


gas and aqueous phases. The results are also presented in Table 6. The results 
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Table 7. Fukui indices of the protonated (HMC) tested inhibitor in aqueous 
medium calculated by DFT/6-311G (d,p). 


ft Í- J+ Í- 
N(1) 0.143 0.002 C(36) 0.029 0.029 
N(2) 0.079 —0.002 C(40) 0.017 0.012 
N(3) 0.064 0.004 C(41) —0.009 0.004 
P(4) 0.017 —0.001 C(44) 0.037 0.032 
P(5) 0.016 0.005 C(48) —0.011 —0.006 
P(6) 0.017 0.004 C(49) 0.019 0.010 
c(7) —0.011 —0.001 C(52) 0.036 0.034 
O(10) 0.028 0.003 C(56) 0.016 0.008 
O(11) 0.023 0.000 C(57) 0.032 0.020 
O(12) 0.026 0.000 O61) 0.026 0.018 
O(13) 0.036 0.002 063) 0.050 0.019 
o(14) 0.009 0.002 O(65) 0.017 0.034 
O15) 0.037 0.003 067) 0.023 0.009 
c(16) 0.013 0.008 O(69) —0.001 0.087 
can) 0.011 0.007 O(71) —0.001 0.138 
C(20) 0.026 0.022 0(73) 0.019 0.083 
c(24) 0.005 0.021 O(75) 0.018 0.116 
C(25) 0.015 0.013 0(77) 0.012 0.060 
C(28) 0.036 0.065 0(79) 0.013 0.050 
CG2) 0.013 —0.002 O81) 0.038 0.009 
CG3) 0.009 0.012 O(83) 0.006 0.067 


indicate that the trend of Ay is also protonated > neutral in both phases, meaning 


that the protonated form is better than the neutral one. 


Finally, Gomez et al [44] proposed the simple charge transfer theory for donation 
and back-donation of charges proposed. In this theory Gomez proved that the inter- 
action between the inhibitors and substrate surface can be probably affected by the 
electronic back-donation process. According to our results given in Table 6, the 
calculated AE,_q exhibits the trend: protonated > neutral in both gas and aqueous 


phases. 


3.6.2. Fukui functions 


The Fukui indices, in corrosion inhibition study, are an important local descriptor to 
investigate the relationship of the inhibitors with mild steel surface in terms of 


structure-activity. Herein, the current research affords an exciting chance to enhance 
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Fig. 8. The optimized geometrical structure of the HMC species. The indicated numbers are used for 


Fukui indices. 


our knowledge of electron donating ability (i.e. nucleophilic fx"), and electron ac- 
cepting ability (i.e. electrophilic f, ) of the inhibitor under probe. The Fukui indices 
of the protonated species of the tested inhibitor in the aqueous phase have been 
calculated based on the Mulliken charges by using Eqs. (14) and (15). The results 
are summarized in Table 7. The indicated numbers of the atoms are given Fig. 8. 
These results agree well the electronic distribution of the FMO orbitals (HOMO 
and LUMO), which were shown in Fig. 9. 


3.7. Monte Carlo simulation 


To get further information and a better understanding of how the HMC molecule 


behave and interact on an iron oxide surface, Monte Carlo simulation were 


EE Oo 
-6.557 x10” - 6.557 x10? eV 


Fig. 9. Optimized geometrical structures, HOMO and LUMO surfaces and electrostatic potential map 
(ESP) of ER (non-protonated) and HMC (protonated). All values of the colored scales are in eV. 
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Fig. 10. Adsorption mode of the HMC on the Fe (110) surface. 


performed on a system with the inhibitor molecule on an iron oxide Fe203 surface, as 
presented in Fig. 10. 


The strength of corrosion inhibitors absorbed on iron oxide surface can be expressed 
by the adsorption energy (Eas). In the present paper, the adsorption energies were 
calculated from the average adsorption energy of the obtained equilibrium configu- 
rations. The highest and lowest obtained F,4, values are -1464.57, -1442.11, kJ/mol 
respectively. A higher negative value of Eg, means that the adsorption of HMC 


molecule is strong, stable and therefore spontaneous adsorption on Fe (110) surface. 


& &- 
H OH 


———=== Electrons transfert n to d 


Electron transfert 1 to d } Donation Chimisorption 
Electron transfert d to z*} Retrodonation 


Fig. 11. Mechanism of HMC adsorption on carbon steel in 3% NaCl solution. 
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As can be seen in Fig. 10, HMC molecule is adsorbed in a parallel orientation to the 
metal surface through the (N, P, O) as active sites [76, 77, 78, 79,]. 


3.8. Inhibition mechanism 


Experimental data, the computational calculations provided a great help in explain- 
ing the mechanism of inhibition of mode of adsorption of inhibitor and identifying 
the active sites. Generally, the charge on the metal, chemical structure (electronic 
density at the donor site) of inhibitor, type of interaction between inhibitor and metal 
influence the adsorption process. The HMC inhibitor has electronegative donor 
atoms such as nitrogen and oxygen atoms and 7-electrons of the aromatic ring. 
These electron-rich donor atoms in addition with 7-electrons of the aromatic ring 
favor adsorption of HMC inhibitor, thus leads to a strong interaction with the steel 
surface. As the observed outcome of the present investigation, HMC under goes 
comprehensive adsorption (involving chemisorption), the possible modes of interac- 
tion of HMC with the carbon steel surface a result of a combination of a three kinds 


of interactions (Fig. 11). 


4. Conclusions 


1. The synthesized Hexa (3-methoxy propan-1,2 diol) cyclotriphosphazene was 
significantly improved the corrosion inhibition efficiency of carbon steel in 


3% NaCl, and their inhibition efficiency increases with increasing concentration. 


2. The results of weight loss (WL) study were consistent with results derived from 
PDP and EIS studies. 


3. The adsorption of HMC on the carbon steel surface from 3 % NaCl obeys a 


Langmuir adsorption isotherm. 


4. The measurement of potentiodynamic polarization demonstrated that inhibitor 
HMC effectively suppressed corrosion process and classified as anodic -type 
inhibitor. 

5. EIS studies revealed the adsorption of inhibitor and are confirmed by increase in 


R, and decrease in Cy values, respectively. 


6. The values of AG a, indicating strongly interaction between inhibitor molecules 


and carbon steel. 


7. DFT computational study shows that the inhibitive strengths of the HMC corre- 
late with their relative reactivity and the presence of highly electronegative 


atoms. 


8. Monte Carlo simulation suggest that the HMC molecules present a parallel 


adsorption model on iron substrate. 
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